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Abstract 

We study the Langevin dynamics of a two-dimensional discrete oscillator chain absorbed on a 
periodic substrate and subjected to an external localized point force. Going beyond the commonly 
used harmonic bead-spring model, we consider a nonlinear Morse interaction between the next- 
nearest-neighbors. We focus interest on the activation of directed motion instigated by thermal 
fluctuations and the localized point force. In this context the local transition states are identified 
and the corresponding activation energies are calculated. As a novel feature it is found that the 
transport of the chain in point force direction is determined by stepwise escapes of a single unit or 
segments of the chain due to the existence of multiple locally stable attractors. The non- vanishing 
net current of the chain is quantitatively assessed by the value of the mobility of the center of 
mass. It turns out that the latter as a function of the ratio of the competing length scales of the 
system, that is the period of the substrate potential and the equilibrium distance between two chain 
units, shows a resonance behavior. More precisely there exist a set of optimal parameter values 
maximizing the mobility. Interestingly, the phenomenon of negative resistance is found, i.e. the 
mobility possesses a minimum at a finite value of the strength of the thermal fiuctuations for a 
given overcritical external driving force. 

PACS numbers: 05.40.-a, 05.60.-k, 05.45.-a, 36.40. Sx 
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I. INTRODUCTION 



Transport phenomena play a fundamental role in many physical systems. For systems that 
evolve in an external potential which possesses metastable states the thermally activated 
escape over potential barriers, as the precondition for transport, is the mostly studied sit- 
uation. The related escape problem, often referred to as the Kramers problem l|, has 
been reviewed e.g. in Due to its ubiquity and simplicity a spatially periodic potentml 
is employed in a number of applications including Josephson tunneling junctions 
phase-locked loops 



jCllljlClJ 



rotation of dipoles, charge-density wave dislocation [8|, diffusion 
of atoms and molecules on crystal surfaces Isl], biophysical processes such as neural activity 
and intracellular transport [lO|, lul, QJfl. Exact expressions for the characteristic quantities 
of motion of one single Brownian particle like the net-current and the diffusion coefficient 
are given in 

In the last decades, the interest in the theory of Brownian motion of interacting particles 
H y. y Q ha^gpwn in .eve.a, fields of science. In this context the stud, of the 



diffusion process [l8| and the mobility of strongly interacting atoms subjected to a 
periodic potential and driven by an external force is a first step towards the understanding 
of solid friction at the atomic level I20I. Recently some studies have considered the 



transport of dimers in a one-dimensional 
spatially uniform dc and ac forces 8, [21I, I22 



p) washboard potential under the impact of 



24 



25| which are applied to all particles. 



A complicated non-monotonous behavior of the mobility depending both on the external 



driving and on the ratio between the competing length scales of the system is found, 
latter plays the role of an internal degree of freedom which generates a ratchet effect 
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32|. 



Particularly in biophysical contexts, the extension to coupled multi-dimensional systems 



e.g. the transport of long and flexible polymers across membranes 



trophoresis 
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33|, 
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35| or DNA elec- 



37l |. has recently attracted considerable interest. Motivated by experimental 
advances in the manipulation and visualization of single polymers using optical and 



magnetic 



39|| tweezers or scanning force microscopy 40!| the external driving acting on the 



system can also be modeled by an external point force which is applied at one single unit 
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In this paper we consider the noise assisted transport of a two-dimensional (2D) discrete 
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oscillator chain confined onto a periodic substrate and subjected to an external localized 
point force. Our theoretical study is related to single-molecule experiments using scanning 
force microscopes. In the commonly used bead-spring model the next-nearest neighbors 
are coupled harmonically. This assumption is valid only for small elongation from the 
equilibrium distance. In order to take the nonlinear character of the coupling into account, 



we introduce an interaction potential of Morse-type which allows bond rupture |44l . |45|. 
Due to the imposed next-nearest neighbors coupling and the fact that the Morse interaction 
potential is rotationally symmetric all configurations with the same distribution of distances 
between two coupled units possess the same amount of energy independent of the distribution 
of the angles between two coupled constituents of the chain. Therefore the bending rigidity 
K and the persistence length Lp = k/T, respectively, are equal to zero. Since the chain is 
stretched during the motion as a result of applied external forces and the substrate potential 
force we can neglect excluded volume interactions. 

In our present study we focus interest on the impact of the ratio between the competing 
length scales of the system and the external driving on the directed motion of the chain. 
The latter is instigated by thermal fiuctuations and quantitatively assessed by the mobility. 
The main question is how both the spatially localized driving and the nonlinear interaction 
potential between the many degrees of freedom infiuence the mobility of the chain? More 
precisely, under which conditions differs the mobility of the considered system strongly from 
the one for ID dimer [3| and does there exit sets of parameters for which the latter 
coincide? 

This paper is organized as follows: In Sec. [Tllwe introduce the model. The transition state 
configurations are identified in Sec. Illli Further, in Sec. [TV]we derive the scaling behavior of 
the activation energy in the limit of weak coupling. The parameter values infiuence on the 
mobility is studied in Sec. [Vl We conclude with a summary and discussion of our results. 



II. THE MODEL 



We study a two-dimensional nonlinear coupled oscillator chain consisting of N particles of 
equal mass m evolving in a 2D substrate under the infiuence of an external dc point force with 
magnitude F. The point force is applied at one single constituent at site uq. The coordinates 
of the n-th unit,n = 1, . . . , A^, in the x-y-plane are given by qn{t) = (a;„(t), yn{t)Y ■ 
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Figure 1: (Color online) Schematic view on the potential landscape with a segment of the chain 
inside. The parameter values are A = 1 and L = 2tt. 

The interaction of the particles v^^ith the substrate is modeled by the on-site potential 

^(a^n) =4 ( 1 - cos ( ) ) , (1) 



2 V V ^ . . 

with periodicity L, i.e. U{x + L) = U{x), and potential height A. Note that the on-site 
potential is translational invariant in the y-direction. The dc point force acting on the unit 
at site no is introduced by the additional potential term 

A segment of the two-dimensional periodic on-site potential U{xn) and the position of the 
chain is depicted in Fig. [U 

Each particle is connected to its two next-nearest-neighbors by nonlinear springs described 
by the Morse potential 

W{rn+i^n) = D [1- exp (-a (r„+i,„ - /q))]^ • (3) 

The Euclidean distance between two units at site n and n + 1 is identified with rn+i,n = 
a/ {xn+1 — XnY + {Vn+i ~ UnY ^hc parameter /q denotes their equilibrium distance. The 
dissociation energy of a bond and the inverse range of the potential are determined by D 
and a, respectively. 

In the overdamped limit the inertia is omitted and the dynamics at finite temperature T 

is described by the Langevin equation (LE) 

1 

IQn = -— Vg-„{W'(r„,„_i) + W{rn+l,n) + U{Xn) -FXn5n,no} + Cn{t) . (4) 
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Figure 2: (Color online) Morse potential W{r) for various inverse interaction range a. The param- 
eter values are D = 10 and Iq = 1. 

Here 7 is the viscous friction coefficient per unit mass. 

For convenient rescaling we introduce suitable space gsc = L/2tt^ energy £^sc = ^, and 
time units tgc = "^7'?sc/-^sc, respectively, and define the dimensionless quantities 

Qsc f-sc -^sc 

F = ^F, D = ^, a=aqsc, I =— . (5) 



Note that the new quantity I determines the ratio between the competing length scales of 
the system Iq/L. Below, we refer to I as the bond length and we omit the overbar in our 
notation. Finally the dimensionless equations of motion read as 

f"n,n—l 



drn,n-l dXn Tn+l 



- - sin (x„) + F 6n,no + at) , (6a) 

^n,n—l 

drn,n-l dyn dv^+l^n 9 Vn 

+ ait) . (6b) 



The stochastic force ^n(t) represents random fiuctuations stemming from the infiuence 
of the environment. The latter is Gaussian white noise with zero mean, {in{t)) = 0, and 
autocorrelation function ( C^'^(t) ^J'^(s) ) = 2T 6ij6x,yS(t — s). 



Throughout this work we impose open boundary condition (OBC), i.e. xi—xo = xn+i — 
xn = and yi — yo = yN+i — Z/Af = 0, and we use an odd number of units N. For the sake 
of symmetry we fix the pulled particle as the one situated at the center of the chain, i.e. 
uq = {N + l)/2. Without loss of generality we set 7, A and m equal to 1 and L = 27i for 
the periodicity. Consequently the scaling parameters equal gsc = Esc = ^sc = 1- 

We remark that the metastable {k even) and unstable {k odd) states of the on-site po- 
tential, dU{xn)/dxn = 0, are located at 

4 = (-l)^arcsin(2F) +fc7r, (7) 

as a result of the external point force. From this follows that in the limit of vanishing 
coupling the pulled oscillator cannot be trapped by the on-site potential given in Eq. ([T]) 
under any circumstances for a tilt larger than the critical value Fcr = 1/2. In order to 
characterize the relative strength of the interaction potential versus the on-site potential, 
we introduce the coupling strength K 

Note that K depends quadratically on the inverse interaction range a and linearly on the 
dissociation energy D. 



III. TRANSITION STATE 



In the following we focus our interest on the escape dynamics of the 2D coupled nonlinear 
oscillator chain from a metastable state over an energy barrier of the corresponding energy 
hypersurface. This progress requires the activation energy -Eact which coincides with the 
height of the energy barrier. We identify the transition state and study the dependence of 
the associated activation energy on the system parameters. 

According to the classical transition state theory [2, 147|, transition states are specials 
points in the 2A^-dimensional phase space. More precisely, a transition state {q^} is a 
hyperbolic fixed point obtained from the stationary system of the deterministic dynamical 
system given by Eq. ([6]) for {i„} = {y„} = 0. The corresponding Jacobian possesses at 
least one real positive eigenvalues which corresponds to the movement along the reaction 
coordinate. All other eigenvalues are negative. We apply a multi-dimensional root finding 
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algorithm using the Newton-Raphson method to calculate the transition state from the 
corresponding stationary system. 

As a result of the considered overdamped dynamics of the system and the reflection of 
the periodicity of the on-site potential U by the 2iV-dimensional phase space instead 
of one global basin of attraction, referred to as the running state, multiple locally stable 
attractors exist [2, llOfl. Therefore the phase space possesses many different transition states 
which can be either linked with one required unique activation energy or different ones. 
Since the escape rate over an energy barrier E'act is assumed to be described by the Van't 
Hoff-Arrhenius law 0, Q|, Tesc oc exp{—Ei^ct/T), we restrict our analytical consideration of 
possible escape processes to scenarios with low amount of activation energy. 

The energy of one configuration is determined by the energy functional V^({q}) 

N 

V{{q})=Yl U{x{) + -Fxn,. (9) 

i=l i=l 

Further, the activation energy i^act is defined as the difference between the energy of the 
transition state configuration {q^} and the energy of the chain at the initial metastable state 

|qmm| 

i^act = V{{c^^}) - V ({q-°}) . (10) 

Below, we discuss the numerically calculated results for the transition state configuration 
(TSC) of the initial escape scenario in which only the pulled particle escapes from its initial 
minima to the next well of the on-site potential in positive x-direction. A sketch of this first 



escape scenario is presented in Fig. 4(a) and is labeled by (1). The numerically calculated 
TSC are depicted in Fig. [3] for various values of the bond length I and for two different 
coupling strengths K. 

At first glance, one can see that for I > x\j, see Eq. ([7l), the escape process is governed by 
an individual escape of the pulled unit independent of the coupling strength. At the TSC 
the riQ— th oscillator is always placed at the top of the energy barrier located at x]j while 
all other units rest at the minimum of the on-site potential t/(x„) without changes of the 
respective bond length /. 

For / < x\j and all values of the coupling strength the pulled particle is placed beyond 
the energy barrier of the on-site potential and its neighbors are elongated in x-direction. 
According to the deterministic part of the equations of motion Eq. (l6bl) . the condition 
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(a) Transition state profiles for weak coupling 
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(b) Transition state profiles for strong coupling 



Figure 3: (Color online) Profile of the transition state belonging to the 1st escape process for 
different values of bond length I and of inverse range parameters a. Only a segment of the chain is 
shown. The dashed line represents the position of the energy barrier xjj of the pulled particle. The 
remaining parameter values are = 99, no = 50,D = 10,F/Fcr = 0.8, a = 0.1 (a), and a = 0.4 
(b). The values of the coupling strength are K = 0.4 (a) and K = 6.4 (b). 

{Vn} = is fulfilled either if the distance between two neighboring particles equals the bond 
length / or the displacement i/n+i — Vn is equal to zero. Since the adopted distance betv\reen 
two neighboring units at the TSC differs from the bond length due to the interaction with 
the substrate U'{xn), the obtained structures resemble the shape of very thin needles, i.e. 
yno±i-yno = 0. The observed refiection symmetry in y-direction, lyno+i-Vml = lym-i-yml, 
results from the choice of hq at the central site of the chain. Note that as a consequence 
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of the imposed open boundary conditions all units of the chain are elongated from their 
starting equilibrium positions at the transition state configuration. 



Comparing the obtained TSC for weak coupling in Fig. 3(a) and the ones for stronger 



coupling in Fig. 3(b) , one recognizes that with increasing coupling the number of elongated 
units involved in the transition state grows. In the case of very weak coupling, the units of 
the chain tend to diffuse via individual steps. Despite that there result large elongations, 
rn,n+i ^ I, during the process the units remain bound to each other. With further increasing 
value of coupling K the oscillators move like a rigid unit with r„,„+i — / and thus display 
organized collective behavior reflected in synchronized escape. 

Such a dependence of the diffusion type on the coupling strength is already experimentally 
known [9|, l49|| and theoretical investigated , l50|| for dimers diffusing on a surface. 

Nevertheless, due to the considered overdamped dynamics of the chain multiple local 
domains of attraction exist for F/ F^r < 1. Hence, the transport of the chain in point force 
direction is determined by stepwise escapes of single units or segments of the chain. These 
stepwise crossings are connected with conflgurational changes of the chain. 



IV. TRUNCATED TRIMER MODEL AND SCALING BEHAVIOR OF THE AC- 
TIVATION ENERGY 

In the following, we study separately the scaling behavior of the activation energy for 
different escape steps in order to assess the time scales of the latter. The character- 
istic time one particle needs to escape from Xq ^ Xq + 2-n (xq is one arbitrary ref- 
erence point) is determined by the first moment of the first escape time distribution 

n n 

Tesc = {t{xQ — > xo + 27r) ) oc exp{Eact/T) [2, 148||. Further the average mean velocity of 
every unit of the chain Vi = 27i/Tesc is determined by Tesc- 



A. First escape scenario 

We start with the first escape scenario. In the latter only the pulled particle escapes from 
its initial minima to the next well of the on-site potential in positive x-direction while 
all other units remain close to their starting position. In order to determine the scaling 
behavior of the activation energy for the first escape scenario in the limit of weak coupling. 
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-C 1, we consider a truncated one-dimensional trimer model. According to the TSC 



presented in Fig. 3(a), the escape process involves only three units for / < x^. All other 
units rest at the minimum of the on-site potential x„ = under maintenance of the bond 
length I and hence do not contribute energy to the activation energy of the escape process. 
The condition of stationarity in y-direction {y„} = is fulfilled by setting = yno±i- 
Since we observed hairpin-like crossing configurations in the simulations for non-vanishing 
but sufficiently point force magnitude F, we consider only reflection symmetry TSC, i.e. 
Xno+i = Xno-1 and i/no+i = Vno-i- FoT I < x^, wc cau assume that the stationary TSC 
is adopted when the pulled unit of the chain is situated close to n and the value of the 
coordinate x„q±i is almost zero. Hence we put 

Xno±l=0 + SXno±l, (lib) 

into the deterministic part of Eq. ( l6al ) and solve the linearized system of equation {in} = 0. 
Despite that we consider the weak coupling limit, <^ 1, further simplifications of the 
system of equations by expanding the interaction potential H^(r„_„_i) up to the 2nd order 
is not possible since the bonds are stretched ,rno,no±i ^ the TSC Fig. [3l Substituting 

A{x) = 2Da (e""("-') - e-^^'^""')) , (12) 
B{x) = 2Da^ (e"-(-') - 2e-='"("-')) , (13) 

we get 



t(i) AB{n)F + 2A{n) 
^no±i- i + 2B{n) ■ ^^^^^ 

For the first escape scenario, the energy of the chain at the initial metastable state is 



almost close to V^^) ({q™''}) = 0.5 (l - Vl - 4F2 - 2F arcsin(2F)) regardless of the value 
of the bond length. Using the expressions given in Eq. (fT4l) . we obtain the following scaling 
behavior of the activation energy for the first escape scenario Eq. ( fTOll 

E^, (xK (^xh-l- sjx\j-l^ (const - F) , (15) 

for K <t^l,l < xIj, and F/F^r < 1. The numerical results for the activation energy for the 



first escape scenario and the analytic one are depicted in Fig. 4(b) and labeled by (1). Both 
results match very well. 
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In the limit I > xlj, the transition state configuration is always given by Xn^^ = xjj 
and Xn^^ = 0,Vn 7^ no, with r^^^-'^i = / regardless of the coupling strength. Then the 
corresponding activation energy attains the saturation value 



El^ = Vl - 4F2 + 2F arcsin(2F) - Fn . (16) 



One recognizes that the curve, label (1), of the activation energy in Fig. 4(b) finally converges 
to the limit value E^l represented by the horizontal dashed line. 

B. Second escape scenario 

In the second escape scenario the no±i— th particles escape from the first well to the next 
one while all other units remain close to their starting position. This case is sketched in 



Fig. 4(a) and labeled by (2). The configuration of the chain at the initial metastable state 
can be described by an effective two particle problem. By using the following ansatz for the 
position of the particles 

Xno = 2tc + 6xno , (17a) 
a^noii =0 + 5xno±i . (17b) 

and solving again the linearized stationary system of equation, we get 

x-(^) - 2vr + 4i?(2vr)F + 4A(2vr)-2F 

^no -^^+ 6i?(27r)-l ' ^^^^^ 

min(2) _ AB{2n)F-2A{2n) 

~ 6i?(27r)-l • ^^^^^ 

The transition state configuration in the 2nd escape scenario is determined by an effective 
three particle problem. For I < xfj — we assume that the stationary configuration is 
adopted when the pulled unit is located close to 27r, its neighbors are located at top of the 
energy barrier, situated at vr, and the no±2— th particles rest at the minimum of the starting 
well. Hence we put 

Xno = 2vr + 6xno , (19a) 
a^noii = TT + 5a;„o±i , (19b) 

Xno±2 = + 5Xno±2 , (I9c) 
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(a) Sketch of probable escape processes 




(b) Activation energy versus bond length 

Figure 4: (Color online) Activation energy Eqx± ei function of the bond length I. The vertical 
dashed line represents x\j and the horizontal dashed line -E^ct given in Eq. (fT6]l . The labels relate 
to the graphs of the activation energy to the corresponding first, second and third escape scenario 
sketched in (a). Note that the saturation value attained by the activation energies -B^);! and 

(3) 

E'act ill the limit of large bond lengths lies above zero represented by the dotted line. The remaining 
parameter values are N = 99, no = 50, D = 10, F/F^r = 0.8, and a = 0.1. 
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into the deterministic part of Eq. ([6]) and solve again the hnearized system of equation under 
the condition {in} = {Vn} = 0. Finally we obtain for the correction terms Sxi the solutions 

,(2) _ (4A(7r)-2F)(l + 2i?(7r)) + 8i?(7r)^F 
^""^^ ~ 4S(7r)2 + 2S(7r)-l ' 

(20b) 



"Oil l + 2B{n) 

I' 



=F- -5xt(2) + 5x2, • (20c) 



In Fig. 4(b), label (2), we depict the graphs showing the dependence of the activation 
energy on the bond length / for the 2nd escape process. The theoretical result matches 
very well with the numerics for / < 27r. In particular it has to be emphasized that the 
position of the minimum of the activation energy as a function of the bond length is very 
well reproduced by the analytic expression given in Eq. (fTSl) and (l20ll . The analytical result 
deviates from the numerics with further increasing value of the bond length since the theory 
does not reproduce that the TSC is always given by Xno^ = x'fj, xll^l_i = it and Xn^"* = 
with Tn^n+i = I for I > xfj — 7T. In addition the activation energy eventually attains the 

(2) 

saturation value -E^ct = 2 in limit of large bond lengths because then the TSC is determined 
by the given above fixed configuration. 



C. Third escape scenario 



At the end of the 2nd escape process the pulled particle and its two neighbors are located 
in the same well of the on-site potential while all other constituents of the chain remain in 



the well behind Fig. 4(a) Therefore the situation for the pulled unit is similar to the initial 
condition in the first escape process and the circumstances for its two neighbors resemble 
the one at the beginning of the second escape scenario. Thereupon, either the no±2— th 
units escape like in the previously mentioned 2nd case or the pulled particle escapes further 



forward. The last case is sketched in Fig. 4(a) and is labeled by (3). The dependence of 



activation energy on the bond length is depicted in Fig. 4(b) The latter is similar to the 



results for the 1st scenario but is evidently shifted to higher energies due to the fact that 
the bonds r„o±i,no±2 are stretched during the escape process which is in contrast to the 1st 
case. Furthermore it is shown that the activation energy reaches the limit value El^l not 
until / x^. 
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To sum up, we have shown that in the limit of weak coupling, K <^ 1, the escape behavior 
of the chain is governed by consecutive individual escapes of single oscillators. In the limit 
of small bond length I <^ x};, the units of the chain escape forward in x-direction by using 
the escape scenario (1) and (2). For larger values of the bond length, the 2nd and 3rd escape 
scenario are utilized by the units in an alternating manner. Notably the vital second escape 
scenario. Sec. IIV Bl requests a higher amount of activation energy compared to the other 
two scenarios and hence governs the time scale of chain transport which is reflected in the 
mobility. 

V. MOBILITY 

Having studied the escape dynamics of the system which is the precondition for transport, 
in the following we focus our interest on the transport properties of a chain of interacting 
nonlinear overdamped Brownian particles conflned onto a periodic substrate. In the case of 
a chain driven by a dc point force oriented along a symmetry axis of a 2D substrate, the 
stationary transport proceeds in force direction, whereas transverse diffusion is not affected 
by the bias. When considering interacting Brownian particles, it is appropriate to study the 
motion of their center of mass (c.o.m.). The position of the c.o.m. in force direction at time 
t is denoted by X{t). Referring to Eq. (iGall . the LE reads as 



With the (5-correlated Gaussian white noise Q{t) = J2i ^ii^)/^- The nonlinear behavior 
of X{t) results from the term Fs = sin (xj) /{2N) in Eq. ( !2Ti l. The latter can be 

interpreted as an effective sliding friction force Fs which results from the interaction of the 
single units with the substrate ?7(x„). The motion of the chain is quantitatively assessed by 
the net velocity of the c.o.m. X{t) 




(21) 



i=l 



Vx = lim 



(22) 



t 



or, equivalently, by the related mobility 



Vx 



(23) 
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In the absence of the on-site potential, the characteristic free net velocity is v^^'^'^ = F/N 
which is associated with the mobility of the chain is /xq = Below, the numerically 

calculated results for the mobility fix are presented in units of the characteristic free mobility 
fio and this is equivalent to express the net velocity in units of the characteristic free net 
velocity v^^""", i.e. fix/ fJ-o = Vx/v^^""". 

The system of coupled LE Eq. ([6]) has been integrated numerically through a second-order 
Heun stochastic solver scheme. Starting from a thermal equilibrated configuration in which 
all constituents of the chain are located near the bottom of one well of f/(x„), the external 
point force has been applied at the no-th unit at to = 0. The stochastic trajectories of the 
constituents have been integrated numerically from to with the time step At = 10^^ ^ tchar 
up to tend = 10^. The characteristic time tchar is given by the relaxation time for the 
overdamped motion of one particle in the biased periodic potential U{xn) in the case of 
vanishing coupling tchar = 2/vT^-4F^. Further average quantities have been obtained as 
ensemble averages over 100 trajectories. 



A. Role of the bond length 

First we study the influence of the internal degrees of freedom, viz. the bond length I, on the 
mobility fix- In Sec. IIV Bl we showed that in the limit of weak coupling the escape rate resc 
and the mean flrst escape time Tesc of the dominating second escape process, respectively, 
exhibits a resonance behavior as a function of the bond length. Due to the fact that escape 
of the units of the chain is necessary for transport of the chain, we expect to observe a 
characteristic dependence of the velocity of the system on the bond length. 

nn 

Several authors [8|,l2j| have shown that the transport of an underdamped ID dimer system 
in a periodic potential strongly depends on its bond length I. In particular, it was found 
that the mobility fix{l) is a reflection symmetric function which attains its maximum value 
at Z ~ (2 A; + l)7r , k E Z, and vari es p eriodically with mod(27r) regardless of the coupling 
strength and force magnitude F sl. l24|. 

Interestingly, it turns out that the mobility of the 2D discrete nonlinear coupled oscil- 
lator chain shows a resonance behavior as a function of the bond length / with one single 
maximum within a period of the substrate potential. This means there exist a set of opti- 
mal parameter values maximizing the mobility and optimizing the transport properties of 
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Figure 5: (Color online) C.o.m. mobility versus the bond length I (a) and the mobility versus the 
mean distance ( r ) (b) for different values of the inverse interaction range a as indicated in the 
plot and with the relation in Eq. jH]) the values of the assigned coupling strength K are 0.4, 3.6, 10, 
and 32.4. The horizontal solid line represents the value of of one single overdamped Brownian 
particle. The remaining parameter values are = 9, no = 5, -D = 10, -F/Fcr = 0.8, and T = 0.1. 

the system, respectively. Further the mobility as a function of system parameter exhibits 
properties which coincide with the results presented in j^, but also shows several new 
phenomena. 

Let us first have a look on the influence of the coupling strength K on the mobility as a 



function of the bond length. In Fig. 5(a) the numerical results for the mobility versus the 



bond length / for various values of a are depicted. We remind that the inverse interaction 
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range a and the coupling strength K are related according to Eq. ([HI). It turns out that 
there exist a finite interval of the bond lengths, l^^^ < I < /^p, in which translocation 
occurs. Further it is recognizable that with increasing coupling strength K the lower limit 
/JP™ approaches from below the position of the unstable state of the on-site potential x\j, 
see Eq. (I7j). This is in compliance with the results presented in Sec. IIV Al since for I > x\j 
the activation energy of the first escape scenario is determined by the saturation value 
Eq. (fT6l) which is independent of the coupling strength K. In contrast in the range I < x]j, 
the number of units participating in the first escape event, and therewith connected the 
activation energy, grows with increasing coupling strength. Hence the mean first escape 
time Tesc oc exp(i?£|/T) increases with the value of K and thus the mobility goes to zero in 
the limit K ^ oo. The upper limit /"p turns out to be independent of the coupling strength. 
Hence the larger the value of K the narrower is the interval of the bond lengths where the 
chain is mobile for a given F and T. 

In addition, it is shown that the mobility grows monotonically with increasing value of 
/ from l^^^ till the bond length is equal to /peak, where the mobility possesses its maximal 
value /i™^''. With upon growing bond length, / > /peak, the mobility decays till / < /^p. It 
turns out that the value of /peak depends on the coupling strength K. In the case of weak 
coupling. K < 1, the posifon eft from . which is in contrast to the results obtained 

for the ID underdamped dimer [8|, l2j|. With increasing coupling strength K the position 
/peak approaches from below vr. 

Further it is demonstrated that the peak height of the mobility /i™^"" increases with the 
coupling strength till K reaches the value Ka- For K > Kcr, the maximum value of the 
mobility remains almost constant upon changing K for a given value of F and T. Note that 
the shape of the mobility graph is asymmetric with respect to /peak which is in contrast to 



the results obtained for the ID underdamped dimer |8l. l24l|. 



Furthermore, the mobility of one sing 
by means of the Stratonovich formula 



e overdamped Brownian particle can be calculated 



161, 



5l| . Comparing the mobility of a monomer, 



represented by the horizontal solid line in Fig. 5(a), with the presented maximum values of 
/i2;(/), one concludes that the response of the extended 2D chain to the external driving is 
less than the one of the monomer. But in the limit of strong coupling both responses nearly 
coincide. This indicates that the chain consisting of nonlinear coupled oscillator synchronizes 
and thus behaves like one heavy Brownian particle with mass A^. 
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Due to the mutual impact of stochastic forces, the point force and the interaction with the 
substrate the bonds of the chain experience dynamical alterations. Therefore the adopted 
distance between two adjacent units mostly differs from the bond length / during the mo- 
tion of the chain. To gain more insight, we discuss the mobility versus the mean dis- 
tance ( r ) between two neighboring units in the long time limit averaged over all sites 



(r) = limt^oo )/(^ Ill Fig. 5(b) the numerical results for the c.o.m. 

mobility fix versus (r) for various coupling strengths are presented. In contrast to the 



results shown in Fig. 5(a), the mobility always reaches the maximal value at (r)peak = tt 
regardless of the value of K. Assuming that the position of a particle Xj+i = Xj + ( r ) can be 
described by means of the averaged distance between two units, one finds that for ( r ) = vr 
the effective sliding friction force Fs, see Eq. ( l2Tll . possesses its minimum value and thus the 
mobility reaches its maximum value. 

Subsequently we study the infiuence of the point force magnitude F on the mobility /j^ 
for given coupling strength and temperature. The results are presented in Fig. [6l It is shown 
that the value of the lower limit value is strongly influenced by the point force magnitude 
F, more precisely the larger the value of F the less is the value of /JP™. Referring to the 
scaling behavior of given in Eq. (fTSl) , the activation energy of the first escape process 



decreases linearly with the value of F for a given bond length /. Hence, the lower limit value 
^low gQgg F/Fcr — > 1. In contrast, it turns out that the upper limit value l^^ is 

not infiuenced by the point force magnitude F . Thus the region where the mobility differs 
significantly from zero strongly depends on F. Similar to the previously discussed results, 
we find that the adopted mean distance at the maximum value of the mobility (r)peak is 
equal to vr. (The corresponding panel is not presented.) Since the deviation of the bond 
length / from ( r ) grows with increasing value of F, one notices that the value of /peak gets 
smaller when enlarging the point force magnitudes F. In addition the corresponding peak 
height increases monotonically with the point force magnitude F. 

Finally, we present the results for the impact of the temperature T on the mobility in 
Fig. [71 For low temperature, ratios of 1/2ti exist where escape and therefore translocation is 
impossible. Upon increasing T, the region around / ~ tt, where the mobility of the chain is 
significantly different from zero, becomes broader. Eventually for high enough T, the chain 
is considerably mobile for all values of I. In addition, it is illustrated that the peak height 
at /peak increases monotonically with T, remaining always smaller than /iq. 
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Figure 6: (Color online) C.o.m. mobility versus the bond length for different values of the point 
force magnitude F . Parameters values are = 9, no = 5, = 10, a = 0.1, and T = 0.1. 
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Figure 7: (Color online) C.o.m mobility versus the bond length I for different values of the temper- 
ature T. The remaining parameter values are N = 9,no = 5, D = W,a = 0.5, and F/F^r = 0.8. 



In Fig. [7] the influence of the bond length on the mobility for a broader interval of values 
of / is depicted. It is shown that ii^ varies periodically with mod(27r). As a novel feature it 
is found that the peak height at /peak decreases upon increasing bond length /. This effect 
is independent of the value of the remaining system parameters. As pointed out in Sec. HVl 
the sequence of occurring escape processes depends on the value of I. For l + n2Ti^ I G [0, 2ti], 
the escape rate r^sc of the flrst n escape processes is determined by r^sc oc exp(— i?^^*/T) 
where El^l is given in Eq. (fTHjl . After the time ti = riT ^ n exp{El^l/T) the subsequent 
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escape processes are similar to the first and second escape scenario respectively for / G [0, 27r]. 
Then the system needs the time t2 ^ n (exp(£^£^(/)/T) +exp(_E^^|(/)/T)) till the next escape 
process occurs similar to the one that has taken place before for / G [0, 27i]. Since the mean 
time till the oscillator escapes, utilizing the 2nd escape scenario, increases with the bond 
length I the peak height at /peak decreases with the latter. 



B. Locked-running transition 

Now we study the impact of the point force magnitude F on the mobility for a fixed value 
of the bond length. In the case of one single overdamped Brownian particle, fi^ and F are 



connected by the Stratonovich formula 
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5l| . A locked-running transition of the system 



occurs if F > with the depinning force F^. The latter is defined in such a way that a 
small change of the external driving strength F = Fd + e,e<t^l, results in a significant 
enhancement of the mobility. 

In Fig. M the dependence of fix on F is presented for different temperatures T. The value 
of the bond length is fixed at / = 3.1 a value for which the chain is considerably mobile 
regardless of the coupling strength K. In general, since one cannot find a preferential direc- 
tion of the random Brownian motion in a spatially periodic potential at thermal equilibrium 



52l | the velocity Vx attains the smallest value for F = independent of the value of T. 
The occurrence of a directed motion Vx would be in contradiction to the second law 
of thermodynamics. Note that for F = the mobility is not defined, see Eq. (l23l) . With 
further increasing value of F the mobility grows monotonically and finally goes to Hq in the 
limit F ^ oo. 

From Fig. [8]we deduce that the value of the depinning force F^ goes to zero with increasing 
temperature T. In particular for T = 0.25, the mobility attains a non-zero value regardless 
of F. For a sufficiently weak point force magnitude, F/F^ < 1, the mobility increases for 
stronger thermal fiuctuations, that is higher temperatures. In this parameter region the 
motion of the chain is mainly instigated by the thermal fluctuations. Additionally, it turns 
out that in the limit of overcritical external driving, F/F^ ^ 1, the slope of the mobility as a 
function of the point force decreases for higher temperatures, i.e. we observe a noise-induced 
suppression of the mobility. 

Since the escape rate and therewith connected the mobility increases with growing value 
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of the temperature according to the Arrhenius-Van't Hoff law, the observed result seems 
to be counter-intuitive. By applying the point force at the riQ-th unit, its effective energy 
landscape changes. The motion of the particles, which are close to the pulled one, proceeds 
preferably in point force direction due to the lower energy barrier. The situation for units 
located nearby the endings of the chain is different. Due to the symmetry of the on-site 
potential the probabilities, respectively, the rates to escape forward or backward are equal. 
Hence it is possible that single units or a segment of the chain escape in direction opposite 
to the one determined by the point force. Due to the fact that the mean time which the 
center of mass needs to move forward a certain distance grows with increasing value of T, 
the effective c.o.m. velocity in force direction becomes lower. 

Finally we study the dependence of the mobility on temperature for certain fixed point 
force magnitudes F. The results are presented in Fig. [9l In the limit of small values of F, 
F/Fcr < 1, the motion of the chain is instigated by thermal fluctuations and thus the c.o.m. 
mobility fix strongly depends on the temperature T. For T = the mobility is equal to zero. 
With further increasing value of the temperature T the mobility fix grows slightly linear. 
In contrast in the limit of a sufficiently strong external driving, F/F^r ^ 1, the motion 
of the chain is purely induced by the external point force F and thus /x^: 7^ for T = 0. 
Upon increasing the temperature the mobility decreases and finally reaches a minimum at 
a finite temperature. By further enhancing the strength of the thermal fiuctuation T, fi^ 
grows slightly linear regardless of the point force magnitude. This phenomenon is called 



negative resistance [53|, l5J, |55||. Comparing the curves for F/Fcr = 1.5 and F/F^r = 2, it 
turns out that the value of the critical temperature depends on the point force magnitude, 
more precisely, the latter is shifted to higher values for stronger point force magnitudes F. 



VI. SUMMARY 



In summary, we have considered the thermal activated motion of an extended two- 
dimensional discrete oscillator chain absorbed on a periodically structured substrate under 
the infiuence of a localized point force. Attention has been paid to the escape dynamics of 
the chain from the metastable states of the substrate potential which is accomplished by 
the adaption of kink-like excitations - also referred to as the transition state. The shape 
of the latter and the corresponding activation energies have been calculated. Due to the 
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Figure 8: (Color online) C.o.m. mobility versus point force magnitude F for different temperatures 
T. The remaining parameters values are = 9, no = 5, D = 10, a = 0.5, and / = 3.1. 




Figure 9: (Color online) C.o.m. mobility versus temperatures T for various point force magnitudes 
F. For the sake of clarity there is a break in the y-axes. Note that the curve for F/F^r = 1, 
represented by the solid line, is a monotonically increasing function of the temperature T. The 
remaining parameter values are = 9, no = 5, L> = 10, a = 0.5, and / = 3.1. 

fact that the bending rigidity is equal to zero, the obtained transition structures resemble 
the shape of very thin needles. As a novel feature we have found that the transport of the 
chain in point force direction is determined by stepwise escapes of a single unit or segments 
of the chain due to the existence of multiple locally stable attractors. In the limit of strong 
coupling, ^ 1, the units exhibit an organized collective behavior and the chain escapes 
like a rigid unit from one domain of attraction to the subsequent one. In contrast for weak 
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coupling, ^ 1, it has been found that the escape behavior of the chain is governed by con- 
secutive individual escape steps of single oscillators. Therefore, we have identified the first 
three possible escape scenarios with the lowest amount of activation energy. Interestingly, 
it has been shown that the vital second escape scenario, which requests a higher amount 
of activation energy compared to the other two scenarios and thus governs the time scale 
it takes for the chain to be transported by one period of the substrate potential, possesses 
a minimum at a certain value of the bond length /. From this follows that for unfavorable 
ratios of the bond length the escape will be highly improbable while other ratios will bear 
good conditions for the chain to escape. The complicated dependence of the activation en- 
ergy, respectively, of the escape rate on the bond length is refiected in a non-monotonous 
behavior of the center of mass mobility. The latter exhibits features which previously have 
been found in models considering one harmonically coupled one-dimensional dimer p, 
viz. the mobility refiects the periodicity of the substrate potential and the latter reaches 
its maximum value if and only if the mean distance between the two coupled units equals 
an odd integer number times the value of the half periodicity of the periodic potential. In 
comparison we have presented several new transport phenomena. For sufficiently weak ex- 
ternal driving, F/Fcr < 1, the mobility shows a resonance behavior as a function of the 
bond length / with one single maximum within each period of the substrate potential whose 
position /peak depends on all system parameters. In general it turned out that the mobility 
attains its maximum value if and only if the mean distance ( r ) between the coupled units 
equals (r)pcak = (2A; + l)7r, A; G Z. The connection between /peak and (r)peak is determined 
by the interplay of different acting forces. Further we have shown that the maximal value of 
the mobility ^"^^^ at /peak decays upon increasing bond length. Lastly, the role of the ther- 
mal fiuctuations and the external driving played for the activated motion of the chain has 
been considered. While it has been found that for sufficiently weak point force magnitude 
the mobility grows monotonically with increasing temperature, for a given overcritical ex- 
ternal driving force the phenomenon of negative resistance has been found, i.e. the mobility 
possesses a minimum at a finite value of the strength of the thermal fiuctuations. 
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